Interplay between function and structure in complex networks 
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We show that abrupt structural transitions can arise in functionally optimal networks, driven by 
small changes in the level of transport congestion. Our results offer an explanation as to why so 
many diverse species of network structure arise in Nature (e.g. fungal systems) under essentially the 
same environmental conditions. Our findings are based on an exactly solvable model system which 
mimics a variety of biological and social networks. We then extend our analysis by introducing 
a novel renormalization scheme involving cost motifs, to describe analytically the average shortest 
path across multiple-ring-and-hub networks. As a consequence, we uncover a 'skin effect' whereby 
the structure of the inner multi-ring core can cease to play any role in terms of determining the 
average shortest path across the network. 
PACS numbers: 87.23. Ge, 05.70.Jk, 64.60.Fr, 89. 75. He 
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I. INTRODUCTION 

There is much interest in the structure of the com- 
plex networks which are observed throughout the nat- 
ural, bi olog ical and social sciences 0j B S S S B Hu 
Ull Il2. I13L n^ . The interplay between structure and 
function in complex networks has become a major re- 
search topic in physics, biology, informatics, and soci- 
ology H M IE E lis 111 El- For example, the very 
same links, nodes and hubs that help create short-cuts 
in space for transport may become congested due to in- 
creased traffic yielding an increase in transit time !6r| . Un- 
fortunately there are very few analytic results available 
concerning network coiigestion and optimal pathways in 
real- world networks 0, 0, S H ES 111 113 ■ 

The physics community, in particular, hopes that a 
certain universality might exist among such networks. 
On the other hand, the biological community knows all 
too well that a wide diversity of structural forms can 
arise under very similar environmental conditions. In 
medicine, cancer tumors found growing in a given organ 
can have very different vascular networks. In plant biol- 
ogy, branching networks of plant roots or aerial shoots 
from different species can co-exist in very similar en- 
vironments, yet look remarkably different in terms of 
their structure. Mycelial fungi jla | provide a particu- 
larly good example, as can be seen in Figs. 1(a) and 
1(b) which show different species of fungi forming net- 
works with varying degrees of lateral connections (anas- 
tomoses). Not only do fungi represent a wide class 
of anastomosing, predominantly planar, transport net- 
works, but they have close parallels in other domains, 
including vascular networks, road and rail transport sys- 
tems, river networks and manufacturing supply chains. 
But given that such biological systems could adapt their 
structures over time in order to optimize their functional 
properties, why do we observe such different structures 
as shown in Figs. 1(a) and 1(b) under essentially the 
same environmental conditions? 

In this paper, we provide exact analytic results for 



the effects of congestion costs in networks with a com- 
bined ring-and-star topology [13 • We thus address the 
question above by showing that quite different network 
structures can indeed share very similar values of the 
functional characteristics relevant to growth. We also 
show that small changes in the level of network conges- 
tion can induce abrupt changes in the optimal network 
structure. In addition to the theoretical interest of such 
phase-like structural transitions, our results suggest that 
a natural diversity of network structures should arise un- 
der essentially the same environmental conditions - as 
is indeed observed for systems such as fungi (see Figs. 
1(a) and (b)). We then extend this analysis by introduc- 
ing a novel renormalization scheme involving cost motifs, 
to describe analytically the average shortest path across 
multiple-ring-and-hub networks. We note that although 
some of the findings of Ref. [ifll appear similar to the 
present ones in terms of the wording of the conclusions, 
the context and structures considered are quite different 
- in addition, the results in the present paper are analytic 
and are obtained in a completely different way. 

As a consequence of the present analysis, we uncover 
an interesting 'skin effect' whereby the structure of the 
inner multi-ring core can cease to play any functional 
role in terms of determining the average shortest path 
across the network. The implication is that any food 
that is found on the perimeter of the network structure, 
can be transported across the structure without having to 
go through the central core - and as a result, the network 
structure in the central core may begin to die out because 
of a lack of nutrient. Interestingly, there is experimental 
evidence that real fungal networks jT^ do indeed exhibit 
such a skin-effect. Other real-world examples in which 
an inner network core ceases to be fed by nutrients being 
supplied from the perimeter, and hence dies out, include 
the vasculature inside the necrotic (i.e. dead) region in a 
growing cancer tumour, and the inner core of a growing 
coral reef. 

Our analytically-solvable model system is inspired by 
the transport properties of real fungi (see Fig. 1(c)). 



A primary functional property of an organism such as 
a fungus, is to distribute nutrients efficiently around its 
network structure in order to survive. Indeed, fungi need 
to transport food (carbon (C), nitrogen (N) and phospho- 
rous (P)) efficiently from a localized source encountered 
on their perimeter across the structure to other parts of 
the organism. In the absence of any transport congestion 
effects, the average shortest path would be through the 
center - however the fungus faces the possibility of 'food 
congestion' in the central region since the mycelial tubes 
carrying the food do not have infinite capacity. Hence the 
organism must somehow 'decide' how many pathways to 
build to the center in order to ensure nutrients get passed 
across the structure in a reasonably short time. In other 
words, the fungus - either in real-time or as a result of 
evolutionary forces ~ chooses a particular connectivity to 
the central hub. But why should different fungi (Figs. 
1(a) and (b)) choose such different solutions under es- 
sentially the same environmental conditions? Which one 
corresponds to the optimal structure? Here we will show 
that, surprisingly, various structurally distinct fungi can 
each be functionally optimal at the same time. 

Figure 1(d) shows our model's ring-and-hub structure. 
Although only mimicking the network geometry of nat- 
urally occurring fungi (Figs. 1(a) and (b)), it is actually 
a very realistic model for current experiments in both 
fungal and slime-mold systems. In particular, experi- 
ments have already been carried out with food-sources 
placed at peripheral nodes for fungi (Fig. 1(e)) and slime- 
mold |16j with the resulting network structures showing 
a wide range of distinct and complex morphologies [13 ■ 
We use the term 'hub' very generally, to describe a cen- 
tral portion of the network where many paths may pass 
but where significant transport delays might arise. Such 
delays represent an effective cost for passing through the 
hub. In practice, this delay may be due to (i) direct con- 
gestion at some central junction, or (ii) transport through 
some central portion which itself has a network structure 
(e.g. the inner ring of the Houston road network in Fig. 
1(f)). We return to this point later on. 



II. THE MODEL 



The Dorogovtsev-Mendes model of a 
Small- World network 




FIG. 1; (a) Typical network for the fungus Resinicium bi- 
colour. (b) Typical network for the fungus Hypholoma fas- 
ctculare. This network has a far denser set of connections 
than (a), yet both are able to thrive in the same environmen- 
tal conditions, (c) Schematic representation of the nutrient 
flows through a mycelial fungus. The food resources (carbon 
(C), nitrogen (N), phosphorous (P)) need to be transported 
efficiently across the entire structure, (d) Our analytically- 
solvable model showing radial connections from peripheral 
nodes to an efi'ective hub. (e) Mycelial fungus Phanerochaete 
velutina after 98 days growing from a central hub-like re- 
source. From day 48, the system is supplied with pairs of 
fresh 4 cm'^ wood resources at 10 day intervals. The resultant 
network has both radial and circumferential connections, as 
in our model (Fig. 1(d)). (f) The man-made road network in 
Houston, showing a complicated inner 'hub' which contains 
an embedded inner ring. 



We begin by introducing the Dorogovtsev-Mendes 
(hereon DM) model [Ijl of a small world network. The 
DM-model consists of a ring-hub structure, and places 
n nodes around a ring, each connected to their nearest 
neighbour with a link of unit length. Links around the 
ring can either be directed in the "directed" model or 
undirected in the "undirected" model. With a proba- 
bility p each node is connected to the central hub by a 
link of length i , and these links are undirected in both 
models. 

We may proceed to solve this model, as in [lj|, by first 



finding the probability P{£, m) that the shortest path be- 
tween any two nodes on the ring is i, given that they are 
separated around the ring by length m. These expres- 
sions can be found explicitly for both directed and undi- 
rected models. Summing over all m for a given i and 
dividing by (n — 1) yields the probability P{t) that the 
shortest path between two randomly selected nodes is of 
length i. The average value for the shortest path across 
the network is then t = X]"=i ^P{^)- For the undirected 



model, the expressions are more cumbersome due to the 
additional possible paths with equal length. However, if 
we define nP{£) = Q{z,p) where p = pn and z = £/n, a 
simple relationship may be found between the undirected 
and directed models in the limit n —^ oo with p —^ 0, that 
is Qundir{z, p) = 2Qdiri'2z , p) . Thus the "directed" and 
"undirected" models only differ in this limit by a factor 
of two: z —^ 2z, with z now running from to 1/2. 



B. The addition of congestion costs 

We generalize the DM model of section Hi Al to include 
a cost, c, for passing through the central hub 13]. This 
cost c is expressed as an additional path-length, however 
it could also be expressed as a time delay or reduction in 
flow-rate for transport and supply-chain problems. We 
then consider a number of cases for the structure of such 
a cost, e.g. a constant cost c where c is independent 
of how many connections the hub already has, i.e. c 
is independent of how 'busy' the hub is; a linear cost c 
where c grows linearly with the number of connections 
to the hub, and hence varies as p = np; or nonlinear 
cost c where c grows according to a number of nonlinear 
cost- functions. 

For a general, non-zero cost c that is independent of £ 
and m, we can write (for a network with directed links): 

p{e,£<c) = -1- (1) 



P{e<m,£>c) = {£-c)p^{l-py''''^ 
P(£ = m,l> c) = 1-p^ J2 («"C)(1 



-P) 



(i-c)-l 



Performing the summation gives: 

P{£ ^ m,£ > c) ^ (l + {£ - c - l)p){l - pY^"-^ 
The shortest path distribution is hence: 



(2) 

(3) 

(4) 



1 

n-l 



y £<c 

+ {n-l-£){£-c)p'^]{l-pY-''-'' \l£>c 



P{£) = <^ 



Using the same analysis for undirected links yields a 
simple relationship between the directed and undirected 
models. Introducing the variable 7 = - with z and p 
as before, we may define nP{£) = Q(z,"f,p) and hence 
find in the limit p — > 0, 71 — » 00 that Qundir{z,^, p) — 
2Qair{2z,2j,p). For a fixed cost, not dependent on net- 
work size or the connectivity, this analysis is straightfor- 
ward. Paths of length I < c are prevented from using 
the central hub, while for ^ > c the distribution P{1) is 
similar to that of Ref. [Tjl . 



III. RESULTS FOR LINEAR AND QUADRATIC 
COST-FUNCTIONS 



For linear costs, dependent on network size and con- 
nectivity, we can show that there exists a minimum value 
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FIG. 2: Our model network showing transport pathways 
through the central hub (connections of length 1/2 denoted 
by thick lines) and around the ring (connections of length 1 
denoted by thin lines). Graph shows average shortest path 
length between any two nodes in a n = 1000 node ring, with 
a cost-per-connection to the hub of fc = 1. There is an opti- 
mal value for the number of connections {p = pn K:i 44) such 
that the average shortest path length i" is a minimum. We 
denote this minimal shortest path length as ^ = ^|min. 
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FIG. 3: Minimal shortest path length ^|miii (i.e. minimum 
value of £) as obtained from Eq. JKJ. (a) Optimal number of 
connections p = pn as a function of the cost-per-connection k 
to the hub. Results are shown for n — 1000 and n = 10000. 
(b) Optimal number of connections p as a function of the 
network size. Results are shown for fc = 2 and fe = 4. 



of the average shortest path £ as a function of the connec- 
tivity to the central hub. We denote this minimal path 
length as ^ = £|min- Such a minimum is in stark contrast 
to the case of zero cost per connection, where the value of 
£ would just decrease monotonically towards one with an 
increasing number of connections to the hub. The aver- 
age shortest path can be calculated from £ = X]"=i ^-P(^) 



from which we obtain 



£{p, n, c 

p(2-2c+2n 



(l-p)"-^(3 + (n-2-c)p) 



1) 



(c-l)(c-7i)pj -3 c(c-l) 



p'^{n — 1) 



2(n- 1) 



(5) 



Figure 2 shows the functional form of i with a cost 
of 1 unit path-length per connection to the hub (i.e. 
c — knp — kp, with fc = 1). The optimal number of con- 
nections in order that ^ is a minimum is approximately 
44 and depends on n. The corresponding minimal short- 
est path £|min is approximately 85. Figure 3(a) shows 
analytic results for the optimal number of connections 
which yield the minimal shortest path ^|min, as a func- 
tion of the cost per connection for a fixed network size. 
Figure 3(b) shows analytic results for the optimal num- 
ber of connections which yield the minimal shortest path 
^Imin, as a function of the network size for a fixed cost 
per connection to the hub. 

To gain some insight into the underlying physics, we 
make some approximations that will allow us to calculate 
the average shortest path analytically for a given cost- 
function that is valid within the approximations. We 
begin by noting that for large n, or more importantly 
large n — c, the first term in Eq. (jsj may be written as 
(1 — p)"~'^ -^ e~P . With the condition that the cost for 
using the hub isn't too high,_the region containing the 
minimum shortest path € = ^jmin will be at sufhciently 
high p to ignore this term, yielding a simplified form for 
the average shortest path: 



p[2 - 2c + 2n - (c - l)(c - n)p\ - 3 



p'^{n — 1) 



c(c — 1) ,„, 



We may then proceed by considering that, for a fixed 
network size and a cost that depends on connectivity, to 
locate the minima we differentiate Eq. © with respect 
to p and set the result equal to zero and obtain 



2 /, , 2 dc 
— (l-c+n) 



p dp dp 



dc , . dc 

c— + (H-n) — + - 



Idc 
2dp 



= 0. (7) 



We substitute into this expression the scaled connectiv- 
ity, p = np, and it then becomes 



^(l-c + n 



n + 



1 



2n\ dc 
p J dp 



(8) 



In the limit of n ^ c and p^ 1, the dominant terms on 
both sides of Eq. © are those in n leaving 



dc 2n 



dp p 



i2 ■ 



(9) 



From this expression we may obtain the location of the 
minimum of the average shortest path for a given cost- 
function for which the approximations are valid. For ex- 
ample, in the case of linear cost c = knp, we find that for 

the optimum number of connections we have p = a / x • 



Using k — 1 and n — 1000, we obtain the value of the op- 
timum number of connections as 44.7, which agrees well 
with the exact value calculated from Eq. Q. Inserting 
the optimum value for p into Eq. © and keeping the 
largest terms we obtain £ « VSfcn, which also agrees well 
with the exact result. 

We now consider quadratic cost-functions, c — kp^. 
This could be a physically relevant cost-function when 
the cost for using the central hub depends on the number 
of connected pairs created, rather than the number of di- 
rect connections made to the hub. Solving for the optimal 
number of connections using Eq. @ gives p - 



responding to a minimum average shortest path length 
£ « \/27kv?. One is also able to consider a cost depen- 
dant on a general exponent c — kp°'. This gives for the 
optimal number of connections 



^"'^ 



2n\ 1+ 



(10) 



The corresponding average shortest path is a more com- 
plicated expression, but it scales with k and n as 



A. 



A:i+°ni+° 



(11) 



This analysis can be adapted to the 'undirected' model 
by using the usual scaling relation between the models 
that was described above. For the case of linear costs 
on an undirected network one gets an optimal number of 
connections at p « ,/?" and a minimum average shortest 

path of £|niin ~ V^kn. 



IV. RESULTS FOR NON-LINEAR 
COST-FUNCTIONS 

We now consider the functional form of £ for non- 
linear cost functions, specifically a cubic cost-function 
and a 'step' cost-function. We show that for these 
non-linear cost-functions, a novel and highly non-trivial 
phase-transition can arise. First we consider the case of 
a general cubic cost-function: 



c(p) = Ap"^ + Bp^ + Cp + D, 



(12) 



where p is the scaled probability, p — pn and A, B,C,D G 
R. In order to demonstrate the richness of the phase 
transition and yet still keep a physically reasonable 
model, we choose the minimum in this cubic function to 
be a stationary point. Hence the cost- function remains a 
monotonically increasing function, but features a regime 
of intermediate connectivities over which congestion costs 
remain essentially flat (like the 'fixed charge' for London's 
congestion zone) . Since we are primarily concerned with 
an optimization problem, we can set the constant D ~ 0. 
Hence 



c{p) = Ap^ - 3Arp^ + 3Ar^p, 



(13) 







FIG. 4: Top: Landscape of the average shortest path length I 
(vertical axis) as a function of the cubic cost-function param- 
eter A and the average number of connections to the central 
hub p. Bottom: The value of p corresponding to a global 
minima in I, as a function of the cubic cost-parameter A. 



where r = ^=^ is the location of the stationary point. 
Substituting into Eq. jSJl yields the shortest path dis- 
tribution for this particular cost-function in terms of the 
parameters A, r, p and n. The result is too cumber- 
some to give explicitly - however we emphasize that it 
is straightforward to obtain, it is exact, and it allows 
various limiting cases to be analyzed analytically. 

Figure 4 (top) shows the value of the average shortest 
path i for varying values of p and A. As can be seen, the 
optimal network structure (i.e. the network whose con- 
nectivity p is such that ^ is a global minimum) changes 
abruptly from a high connectivity structure to a low con- 
nectivity one, as the cost-function parameter A increases. 
Figure 4 (bottom) shows that this transition resembles 
a first-order phase transition. At the transition point 
A = Acrit, both the high and low connectivity structures 
are optimal. Hence there are two structurally inequiva- 
lent networks having identical (and optimal) functional 
properties. As we move below or above the transition 
point (i.e. A < Ac^a ov A> Acnt respectively) the high 
or low connectivity structure becomes increasingly supe- 
rior. 

Using the same approximations as those in Sec. IIIII we 



may estimate the approximate value of the cubic function 
parameter, A, that will lead to two minima in the average 
shortest path distribution. We proceed by solving Eq. 
(jSJ with our cubic function as the cost-function: 

dc 
p2— -2n^ 3A(p2 - 2rp + r'^)p^ - 2n = . (14) 
dp 

The solutions to this equation are then stationary points 
in €, and at least three stationary points are required for 
the distribution to have multiple minima. We thus have 
p = and p = (3 ± \/7)r. Inserting the central value, 
(3 — V7)r, into Eq. H14|l gives an approximate lower 
bound for A: 



A„ 



n 



(15) 



Although this calculation does not give us the value of 
Acrit, it is expected (and results confirm such a conjec- 
ture) to be close to Amin- From this analysis, we can also 
see that both the location of the minima and the distance 
between them is governed by the cubic parameter r. 

We have checked that similar structural transitions can 
arise for higher-order nonlinear cost-functions. In partic- 
ular we demonstrate here the extreme case of a 'step' 
function, where the cost is fixed until the connectivity 
to the central hub portion reaches a particular threshold 
value. As an illustration, we consider the particular case: 



c(p,ro) = 50 ^Sgn(p 



50 



^ro 



50 



(16) 



where Sgn(a;) = — 1, 0, 1 depending on whether x is nega- 
tive, zero, or positive respectively, and rg determines the 
frequency of the jump in the cost. Figure 5 (top) shows 
the average shortest path I for this step cost-function 
(Fig. 5 (bottom)) as p and rg are varied. A multitude of 
structurally-distinct yet optimal network configurations 
emerge. As tq decreases, the step-size in the cost-function 
decreases and the cost-function itself begins to take on a 
linear form - accordingly, the behavior of £ tends towards 
that of a linear cost model with a single identifiable min- 
imum. Most importantly, we can see that once again a 
gradual change in the cost parameter leads to an abrupt 
change in the structure of the optimal (i.e. minimum t) 
network. 



V. THE RING-HUB STRUCTURE AS A 
NETWORK MOTIF 

We have allowed our ring-and-hub networks to seek 
optimality by modifying their radial connectivity while 
maintaining a single ring. Relaxing this constraint to 
allow for transitions to multiple-ring structures yields a 
number of related findings. In particular, allowing both 
the radial connectivity and the number of rings to change 
yields abrupt transitions between optimal networks with 
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FIG. 5: Top: Landscape of the average shortest path length 
£ (vertical axis) as a function of the 'step' cost-function pa- 
rameter ro and the connectivity p. Bottom: The 'step' cost- 
function as a function of the step-frequency parameter ro and 
p. As ro decreases, the cost-function becomes increasingly 
linear. 



different radial connectivities and different numbers of 
rings. One could, for example, consider this to be a model 
of a complicated fungal structure (Figs. 1(b), 1(e)) or of 
the interactions between hierarchies within an organiza- 
tion, such as in Fig. 6. 

To analyze analytically such multiple-ring structures 
we introduce the following renormalization scheme. Con- 
sider the two-ring- and- hub network in Fig. 1(f). For 
paths which pass near the center, there is a contribution 
to the path-length resulting from the fact that the inner 
ring has a network structure which needs to be traversed. 
Hence the inner-ring-plus-hub portion acts as a renormal- 
ized hub for the outer ring. In short, the ring-plus-hub 
of Eq. (1) can be treated as a 'cost motif for solving 
multiple- ring- and- hub problems, by allowing us to write 
a recurrence relation which relates the average shortest 




FIG. 6: Schematic description of hierarchies in a human or- 
ganization, institution or company. As shown, this diagram 
can be re-drawn as a multiple-ring-and-hub structure. Similar 
networks are likely to exist in a range of social systems. 



path in a network with i + 1 rings to that for i rings: 

it+i{Pt+i,ni+i,c) = 
(1 - p,+i)".+i-''.(p.."..'=) (^3 + (n,+i - 2 - lipi, m, c))pi+i) 



+ 



2pj+in - ii(pi,nt,c) + ni+i) 
vi+Ani+i -1) 



Pi+i((^i(Pi,ni,c) - l)(^i(pi,n,;,c) - ni+i)pi+i ) - 3 



pi+An-i+i -1) 
£i{pi,ni,c){li{pi,nt,c) - 1) 
2(n,+i - 1) 

where i > and £o — c with c being a general cost for 
the inner-most hub. The case i = is identical to Eq. 
(1). As before, pi+i represents the probability of a link 
between rings i + l and i and rii+i is the number of nodes 
in ring i + l. 

We investigate the properties of our renormalized TV- 
ring network by placing a number of constraints on the 
parameters and observing the average shortest path be- 
haviour, so we may once again determine regimes of func- 
tionally optimal network configurations. We begin by 
increasing the number of rings to iV = 2, with the con- 
straint that the number of nodes on each ring is fixed. We 
find that the configuration that yields the minimum av- 
erage shortest path length has all the ring probabilities. 
Pi, equal. Figure 7 demonstrates the average shortest 
path distribution for such a case. Figure 7 also demon- 



strates the accuracy of our analytic renormalized result, 
as compared to a full-scale numerical calculation for £. 
If we allow the number of nodes on the two rings to be 
different, we find that the configuration that optimizes 
the shortest path favours a larger number of connections 
on the ring with the most nodes. Returning to the orig- 
inal configuration, N = 2 with an equivalent number of 
nodes on each ring, we consider the effect of a cost on 
the central hub of the inside ring. We find that a greater 
number of connections on the ring without costs opti- 
mizes the network, as one might expect. 
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FIG. 7: Average shortest path £ for a network with N = 2 



rings, each of size n 
P2 =p. 



10 and with equal connectivity pi 



The addition of further rings, such that N > 2, leads 
to some interesting results. For a network with an equiv- 
alent number of nodes on each ring the optimal configu- 
ration remains such that all the ring probabilities, pi, are 
equal. However, for N ^ 2, we begin to see a deviation: 
connections should be moved to the outer rings of the 
network (those furthest away from the hub) in increasing 
numbers out to the edge of the network to obtain the 
minimum shortest path. We thus consider the properties 
of a network with varying N; equal ring probabilities pf, 
inner ring costs, cq, that vary from — > oo; both fixed 
numbers of nodes on each ring, ni = n — constant, and 
varying numbers of nodes on each ring such that Ui cc i. 
Figure 8 shows the average shortest path distribution for 
several such cases. Interestingly, for large N the system 
becomes indifferent to the cost of the central hub, cq, 
and all the distributions converge, in both the case of 
fixed and varying numbers of nodes on each individual 
ring. This is suggestive of an effective 'skin effect', as the 
center of the network becomes effectively disconnected 
after the addition of a large number of rings. The re- 
dundance of the central portion of the network offers an 
explanation for an earlier finding: that as we approach 
N ^ 2 the optimal configuration ceases to be such that 
all innermost ring probabilities, pi, are equal. We find 
that to optimize the network we need to move the con- 
nections into the skin, more than likely as a result of the 



detachment of the center of the network from the whole. 
By comparison of the shortest path values of our mul- 
tiple ring networks, we have found a further important 
result. We find that there are optimal network structures 
with different numbers of rings and radial connectivities, 
yet which have the same average shortest path length 
across them |l7l |. Hence, as before, optimal network 
structures exist which are structurally very different, yet 
functionally equivalent. Figure 9 shows an explicit exam- 
ple of two such functionally equivalent, optimal networks. 
It is remarkable that these images are so similar to the 
real fungi shown in Figs. 1(a) and (b). 



VI. CONCLUSION 

In summary, we have analyzed the interplay between 
the structure and function within a class of biologically- 
motivated networks, and have uncovered a novel struc- 
tural phase transition. Depending on the system of inter- 
est (e.g. fungus, or road networks) these transitions be- 
tween inequivalent structures might be explored in real- 
time by adaptive re- wiring, or over successive generations 
through evolutionary forces or 'experience'. Through the 
use of an approximation, we treated the original network 
as a cost-motif for a larger network and considered the 
circumstances under which such networks obtained their 
optimal functional structure. The equivalence in func- 
tion, defined by the networks transport properties, be- 
tween various topologically distinct structures may pro- 
vide insight into the existence of such disparate structure 
in real fungi. An important further implication of this 
work is that in addition to searching for a universality in 
terms of network structure, one might fruitfully consider 
seeking universality in terms of network function. 



Acknowledgements 

We kindly acknowledge L. Boddy, J. Wells, M. Har- 
ris and G. Tordoff (University of Cardiff) for the fun- 
gal images in Figs. 1(a), (b) and (e). N.J. and M.F. 
acknowledge the support of EU through the MMCOM- 
NET project, and N.J. acknowledges the support of the 
EPSRC through the Life Sciences Program. 



g 

1"^ 



10* 



: (a) 


' ' ' ' 


; 


— = 10"= 

— c =61215 

c =0 


- 


r\^_ 




-: 


r ^ 




N 



— 


c 


= 10« 


— 


c 


= 10= 


___ 


c 


= 2x10'' 


— 


c 


= 


10^ 


-- 


,11,1 , 
10^ 



10* 



AT 



FIG. 8: Average shortest path i" as a function of the number of 
rings TV that make up the network. In (a) all rings are of size 
n = 10® and in (b) the rings increase according to rii — lO^i 
where i is the ring index. In both cases, the probability of 
connectivity of a node is a constant and equal to pi — 0.001. 
The innermost ring has a hub with a constant cost co . Here cq 
ranges from to effectively infinite. In all cases the limiting 
value for £ is the same, demonstrating the 'skin effect' and 
hence the effective disconnection of the inner rings. 




FIG. 9: Two structurally inequivalent networks, which are 
functionally equivalent in the sense that they have the same 
average shortest path across them. The average shortest path 
across the structure is the same for both networks, and both 
networks are themselves optimal (i.e. minimum £). Nutrients 
found on the perimeter of each structure should therefore take 
the same average time to cross it - this in turn implies that 
both structures would co-exist in the same favourable envi- 
ronmental conditions, (a) A single-ring-and-hub network with 
a linear cost-function, (b) A two-ring-and-hub configuration. 
The inner ring-hub structure has the same cost-function as in 
(a). The similarity to the real fungi in Figs. 1(a) and (b) is 
striking. 
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